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FIGURES 


Variation  of  Re  cm  and  Im  cm  with  R for  the  first  two 
discrete  temporal  eigenvalues  of  the  model  equation 
when  


Variation  of  Re  <*m  and  Im  with  R corresponding  to 

waves  propagating  downstream  for  the  first  two  distrete 

spatial  eigenvalues  of  the  model  equation  when  <*>  = 1 

Variation  of  Re  and  Im  Qfm  with  R corresponding  to 

waves  propagating  upstream  for  the  first  two  discrete 

spatial  eigenvalues  of  the  model  equation  when  u>  = 1 

Variation  of  Re  a and  Im  or  with  R corresponding  to  waves 
propagating  downstream  for  the  first  two  discrete  spatial 
eigenvalues  of  the  Or  r - Somme  rf  eld  e /nation  when  = 0.  0397  . . 
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1.  INTRODUCTION 

The  recent  studies  by  Jordinson*  and  Mack2  of  the  temporal  eigenvalues 
of  the  Blasius  boundary  layer  have  provided  new  insight  into  their  properties 
and  raised  a number  of  questions  about  their  role  in  initial  value  and  other 
problems.  Mack  found  that  for  a given  Reynolds  number  and  wave  number 
there  exists  a finite  number  of  discrete  eigenvalues;  as  the  Reynolds  num- 
ber is  increased,  additional  eigenvalues  spring  from  the  continuous  spectrum. 
Mack  discusses  the  properties  of  the  continuous  spectrum,  which  have  also 
been  considered  by  Grosch  and  Salwen  (unpublished).  These  studies  show 
that  the  eigenfunctions  of  the  continuous  spectrum  vary  sinusoidally  as  y goes 
to  infinity.  Thus  these  functions  are  bounded  but  are  not  zero  in  this  limit. 
Our  purpose  in  this  paper  is  to  offer  an  interpretation  of  the  continuous 
spectrum  and  an  explanation  of  the  properties  of  the  discrete  spectrum.  We 
shall  also  make  a parallel  study  of  the  spatial  eigenvalues. 

2.  A MODEL  EQUATION 
Consider  the  differential  equation 

with  0=0  when  y=0,  0 — 0 as  y — oo.  Here  R^  is  a (constant)  Reynolds 
number  and  U(y)  a given  velocity  distribution  taken  to  be  U = 0 when  y < 1, 

U = 1 when  y > 1 . 
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We  are  primarily  concerned  in  this  paper  with  the  spectral  properties 
of  the  Or r - Somme rfeld  equation  with  a Blasius  velocity  profile.  The  model 
equation,  Eq.  (1),  is  introduced  because  its  spectral  properties  are  generally 
similar  to  those  of  the  Or r -Sommerfeld  equation  and  as  a result  we  are 
better  able  to  understand  the  latter.  A particular  advantage  of  the  model 
equation  is  that  the  eigenvalue  problem  can  be  reduced  to  a simple  transcen- 
dental equation.  A potential  problem  with  the  model  equation  arises  from  the 

use  of  a discontinuous  U,  a profile  which  might  introduce  spurious  modes 

2 

into  the  spectrum.  Mack  has  found  that  the  discrete  spectrum  of  the  Orr- 
Sommerfeld  equation  changes  significantly  when  the  velocity  profiles  are 
changed  from  analytic  to  nonanalytic  functions.  As  will  be  seen,  however, 
the  model  equation  with  the  nonanalytic  profile  has  asymptotic  character  and 
spectral  properties  similar  to  the  Or r -Sommerfeld  problem. 

For  temporal  eigenvalues  we  write 

<f>  - $(y)  exp(iax  - iacmt)  (2) 


where  a is  real,  c is  the  eigenvalue  to  be  found  and  $,  are  taken  to  be 
m 

continuous  at  y = 1.  We  shall  show  that  a discrete  spectrum  exists  with 
eigenfunctions  which  are  zero  at  infinity  as  well  as  a continuous  spectrum 
with  bounded  eigenfunctions  at  infinity. 

Consider  first  the  discrete  eigenvalue  problem  which  may  be  reduced  to 
the  solution  of 

e2  = - 6i  coth  0j  (3) 
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where 


G2  = a2  - i ac  R,  02  = a2  + iaR(l-c  ) 
1 m dd  m 


(4) 


and  ReG^  > 0*  When  R » 1 it  is  easy  to  show  that 

1 \T2 
172 


-i  , 2 . 2 2,  , n2TT2  \T2  /4  , .4  , 
cm  ~ aR  + n u ) + T7 T (1  + i)  + • • • 


(5: 


(oR) 


where  n is  an  integer.  Thus  in  the  limit  R-»oo  there  are  an  infinite  number  of 

discrete  eigenvalues.  It  may  be  shown  that  there  are  no  discrete  eigenvalues 

for  sufficiently  small  values  of  R and  that  Re  c * 0,  Im  c ^ 0.  Denoting 
’ mm 

th 

the  n eigenvalue  by  cmn>  we  nc>w  consider  its  variation  with  R,  assuming 

2 

it  to  be  continuous.  By  analogy  with  Mack's  numerical  work,  we  expect  the 

discrete  eigenvalues  to  terminate  in  the  continuous  spectrum  and  therefore 

look  for  a value  R of  R such  that  Re  c is  unity  (as  will  be  seen  subse- 
n mn 

quently,  Re  c^  = 1 is  in  the  continuous  spectrum).  At  this  value  of  R we  set 


mn 


= 1 - , H-2  = (v-l)a2,  Gj  = a-ib,  02  = «+ifx  (6) 

n 


where  v,  a,  b,  e are  all  real  and  positive.  The  quantity  € is  arbitrarily 
small  and  can  be  neglected  (we  consider  points  arbitrarily  close  to  the  con- 
tinuous spectrum).  With  further  manipulation,  we  obtain 


b2  = ^2  + a2,  ab  = aR  /2 

n 


(7) 


It  may  readily  be  shown  that 


b = a cosh(2a),  a = -bcos(2b),  sin(2b)  = -tanh(2a),  p = a sinh(2a);  (8) 
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again  there  are  an  infinite  number  of  solutions  of  which  the  smallest  has 
a = 0.815,  b = 2.163,  aRj  = 3.526,  cm  = 1 - i (0.  284a^  + 1.138) 
and,  when  n is  a large  integer 

■ 3 2 

11  1Q  Rn 

b *(n  - 4)Tr’  a*log[(n-TM/2,  «Rn  = 2ab,  * 1 - 


(9) 


Further,  it  may  be  shown,  and  is  in  fact  implicit  in  Eq.  (8),  that  there  are 

no  solutions  in  which  p is  either  negative  or  imaginary.  In  Fig.  1 we  display 

the  variation  of  Re  c and  Im  c for  two  of  the  principal  eigenvalues. 

mn  mn  r r 

The  inference  is  clear.  When  R takes  an  infinite  value  there  is  also  an 

infinite  number  of  discrete  eigenvalues.  At  finite  values  of  R there  is  only 

a finite  number,  and  for  each  there  is  a critical  value  R of  R at  which 

n 

Re  c reaches  unity;  at  smaller  values  of  R this  eigenvalue  is  no  longer 
mn 

physically  meaningful.  An  analytic  continuation  of  the  solution  of  Eq.  (3) 
may  be  found  when  R < R , corresponding  to  Re  c >1,  but  then  Re  8,  < 0 
and  therefore  we  have  a problem  in  which  it  is  desirable  to  remove  the  expo- 
nentially decaying  solution  as  y— *00.  It  is  natural  to  ask  whether  there  are 
any  other  solutions  in  which  Re  c >1  and  Re  6.,  < 0;  the  answer  appears  to 
be  no.  Certainly  if  such  solutions  intersect  the  continuous  spectrum,  then 

as  Re  c — 1 + we  have  8-,  = ip  with  p < 0,  and  we  know  that  no  such  solutions 
m c. 

exist. 
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Figure  1 


when  a = 1.  They  branch  off  the  continuous  spectrum 
when  R =3.53  and  12.  1,  respectively. 
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If  we  admit  solutions  which  remain  bounded  as  y-»oo  but  do  not  neces- 


sarily tend  to  zero,  the  continuous  spectrum  exists  for  all  R. 
tions  are 

c = 1 - iaN/R,  N > 1,  = - ior  - Qf?'(N-l),  0?  = 

m \ c. 

$>  = 0^ sinh(e^y)f  y < 

1 /2 

$ = 0^sinh  0^cos  [er(y-l)(N-l)  ] 

+ 02  cosh  OjSin  [cy(y-  1 ) (N - 1 ) 1 7 2 ] / [cr (N - 1 ) 1 7 2 ] , 


These  solu- 


io-  ( N-  1) 1 /2 

(10) 


y > 1 (11) 


However,  another  interpretation  can  be  placed  on  these  solutions  which  is 
more  satisfactory  from  a physical  standpoint.  Suppose  that  at  time  t = 0 

<f>  = f(y)e^QX  (12) 

where  f(y)  is  a given  function  of  y vanishing  at  y = 0 and  in  the  limit  as  y — *-oo. 
How  does  $ behave  at  large  values  of  t,  where 

<t>  = ®(y,  t)eiaX,  - R |^  = (a2  + iaUR  )$>  (13) 

dy 

It  then  consists  of  the  sum  of  a set  of  eigenfunctions  satisfying  null  boundary 
conditions  at  y - 0,oo  and,  if  there  are  any  discrete  eigensolutions , i.e.,  if 
R > R these  must  also  be  included.  The  coefficients  of  such  eigenfunctions 
may  readily  be  determined,  for  example  by  Laplace  transform  analysis,  in 
terms  of  f(y).  In  addition  there  are  always  eigenfunctions  derived  from  the 
fact  that,  for  y > 1,  the  governing  equation  is  of  the  heat-conduction  type. 


Mnmn 


Therefore,  we  may  expect  that  at  large  values  of  t the  solution  of  Eq.  (13)  in 
y > 1 will  include  an  expression  of  the  form 


. 2 D - 1 . „ 1/2  -1/2. 

? = exp(-  a tR  - iat)  J t F^(yR  t ) 


due  to  this  cause.  Here  satisfies  Hermite's  equation,  e.g., 

2 

F = exp(-y  R/4t)  and  F^-0  as  y— oo.  At  y = 1 this  solution  and  its  first 
derivative  with  respect  to  y must  be  continuous  with  a solution  of  Eq.  (13)  in 
0 < y < 1.  This  solution  may  be  expressed  in  the  form 

2,„  - 1 . , . -P  , , v /«  C V 


$ = exp(-  a tR  - iat)N  ' t 13  f (y) 

o r j P 


where  f is  a complex  exponential  function  satisfying  f (0)  = 0.  By  matching 
P P 

Eqs.  (14)  and  (15),  across  y = 1 we  establish  immediately  that  F (0)  = 0 for 
them  conjointly  to  describe  tbe  eigenfunctions  as  t— -oo.  This  condition  in 
turn  implies  that  m must  be  a positive  integer  and  p = m + 1/2. 

Thus,  we  see  that  at  a finite  value  of  R the  eigenfunctions  consist 

mainly  of  an  infinite  set  of  functions  which  decay  exponentially  in  time,  decay 

1 /2 

exponentially  with  y outside  a region  of  thickness  (t/R)  , and  have  a 

boundary-layer  structure  in  0 < y <1.  Insofar  as  the  effect  on  the  external 
region  y > 1 is  concerned,  these  are  the  most  significant  because  their 
thickness  grows  with  t.  The  phase  velocity  of  these  solutions  is  given  by 


c = 1 - ia/R.  Now  each  of  the  functions  F can  be  expressed  in  the  form 
m m 


C 8m(0)exP  - i0y  - Q^tR  d0 


and  so  may  be  regarded  as  a sum  of  the  eigenfunctions  of  the  continuous 

spectrum  (Friedman  has  made  such  an  identification  for  the  continuous 

spectrum  of  the  membrane  equation).  In  addition  to  the  parabolic  eigen- 

solutions  there  is,  for  R > R,,  another  discrete  set  of  the  form  of  Eq.  (2) 

in  which  0 < Re  c <1.  As  R increases,  this  set  is  augmented  at  discrete 
m 

values  of  R by  new  eigensolutions , which  may  be  regarded  as  springing  from 
the  continuous  spectrum  whenever  Eq.  (3)  is  satisfied,  so  that  at  these  points 


Eq.  (11)  reduces  to 


- 0(  sinh  0^  exp  - ia(y-l)  J N-  1 in  y > 1 


3.  TEMPORAL  SPECTRUM  OF  THE  BLASIUS 
BOUNDARY  LAYER 

The  Orr -Sommerfeld  problem  may  be  formulated  as  follows.  Consider 
the  equation  for  the  small  perturbations  of  the  straamfunction  for  a steady 
shear  flow  defined  by  a velocity  U(y)  in  the  x direction.  Disregarding  the 
fact  that  the  steady  flow  does  not  satisfy  the  Navier -Stokes  equations,  we 


have 


^-(vz4d  + U^(vZ*)  - U"(y)^  - R'V* 


where  R^  is  the  Reynolds  number.  For  the  Blasius  boundary  layer  U(0)  = 0 

and  U(oo)  = 1.  Further,  xx^R  * measures  distance  along  the  plate  from  the 

center  of  the  disturbance,  itself  a distance  xq  from  the  leading  edge  of  the 

plate;  yx  R * measures  distance  normal  to  the  plate;  tx^R  * measures  time; 
2 

and  R is  based  upon  x . 

The  boundary  conditions  satisfied  by  may  be  taken  as 
vb 

4»  = r — =0  at  y = 0 and  ^—0  as  y— oo  (19) 

oy 

for  all  x and  positive  t.  The  assumptions  leading  to  Eq.  (18)  are  valid  when 
R » 1,  but  here  we  shall  consider  the  properties  of  the  solufion  of  Eqs.  (I8i 
and  (19)  for  all  R. 

The  standard  method  for  finding  the  temporal  eigenvalues  of  is  to 

write 

+ = ei0,x-cl)t(yl  (201 

with  a being  real,  whereupon  'k  satisfies  the  Or r -Sommerfeld  equation  with 
c as  eigenvalue.  This  is  a fourth-order  equation  and,  as  y— *00,  has  a richer 
asymptotic  structure  than  our  model  equation  which  is  only  second  order.  In 
addition  to  the  term  e [cf.  Eq.  (4)],  ^ can  include  a term  a e a^.  The 

model  equation,  if  y > 1,  is  actually  the  same  as  that  for  ^ if  y » 1,  so  it  is 
not  surprising  that  in  spite  of  its  relative  simplicity,  practically  all  the 
general  properties  of  <$  and  c have  their  parallel  in  those  of  'if  and  c as 


elucidated  by  Mack.  First  of  all  there  is  a continuous  spectrum  given  by 

c = 1 - (p2  + a2),  P real,  (21) 


for  which  ik  is  merely  bounded  as  y— co. 

This  spectrum  may  also  be  interpreted  in  terms  of  the  discrete  spec- 
trum of  the  heat  conduction  equation.  We  write 


i,  ia(x-t)  - a t/R  . . 

4*  = e (y,  t) 


(22) 


and  set 


oo 


1/2 


^2  ^n(y)t  n i£^‘' 


n = 0 


(23) 


where  m,  p are  real,  p - 0 and  p > p ..  Then  'I'  satisfies 
“n  o n n - 1 o 

- a2*"  + iaRU"(y)*  = 0 (24) 


and  has  a solution  in  which  ^(0)  = ^(0)  = 0 and 


•q/  ~Ae°^  + By  + C 
o o o o 


as  y— -ao 


(25) 


There  is  no  numerical  evidence  that  the  constant  Bq  is  ever  zero.  When 

1/2 


y ~ t and  t is  large,  we  write 


* = t' 


oo 


n = 0 


2~ 


$ (y)t  ?n  ' + F (yR1/2t‘1/2)t  9,1 

n n 


(26) 


where  a > a , , *n  - a w = 0 and  F is  defined  in  Eq.  (14).  A match 
^n  ^n-1  o o o 


-14- 


M 


between  Eqs.  (25)  and  (26)  is  only  possible  if  ^ = A^e  = 0, 

F'(0)  = R , and  q =0.  Since  F^  — 0 as  y/t*^— oo,  it  follows  that  m is 

an  integer.  Thus,  as  with  the  model  equation  there  is  a discrete  spectrum 
of  the  heat-conduction  type,  each  member  of  which  may  be  interpreted  as  an 
integral  over  the  continuous  spectrum. 

The  eigenfunctions  [Eqs.  (23),  (26)]  should,  from  a s trie t mathematical 

point  of  view,  be  regarded  as  the  central  set  of  eigenfunctions.  In  addition, 

for  given  a,  there  is  a critical  value  Rq  of  R beyond  which  a second,  discrete 

set  springs  from  the  continuous  spectrum.  This  behavior  is  similar  to  that 

2 

of  the  model  equation.  The  results  of  the  two  searches  performed  by  Mack 

and  the  authors  suggests  that  at  the  intersection  points  of  the  continuous  and 

discrete  spectrum  'I',  which  oscillates  finitely  when  y » 1,  has  a form  similar 

to  Eq.  (17),  especially  in  regard  to  the  sign  of  the  argument  of  the  exponential. 

No  solutions  were  found  in  which  this  sign  was  positive.  This  suggests  that 

there  is  no  branching  from  the  continuous  spectrum  leading  to  eigenvalues  in 

2 

which  Re  c > 1.  The  value  of  R was  shown  by  Mack  to  be  15  when  0=0.  179 

o 1 

and  we  found  it  to  be  25  when  a = 0.115.  As  R increases,  more  members  of 
this  second  set  appear  and  the  total  number  tends  to  infinity  with  R.  The 
properties  of  these  eigenvalues  appear,  with  one  exception,  to  be  consistent 

with  those  of  the  model  equation;  i.e.,  with  Re  c— *0  + , Im  c — 0-  as  R— oo.  The 

2 2 
exception  is  the  inviscid  eigenvalue,  which  Mack  has  labelled  3 when  R ~ 10  , 

as  shown  in  Fig.  3 of  his  paper.  For  this  one  c tends  to  a finite  limit.  These  two  types 

, of  limiting  behavior  are  consistent  with  the  two  classes  of  asymptotic  solutions 


Ik  „ 


described  by  Lin.  A farther  special  feature  of  these  eigenvalues  distinguish- 


ing them  from  c^,  and  of  crucial  importance,  is  that  their  real  and  imaginary 


parts  are  not  necessarily  monotonic;  it  is  possible  for  Im  c > 0. 


4.  THE  CONTINUOUS  SPECTRUM  OF  SPATIAL 
EIGENVALUES 


The  spatial  eigenvalues,  obtained  by  requiring  R and  w = ac  to  be  real 


quantities  and  looking  for  (possibly  complex)  values  of  a which  permit  non- 


trivial solutions  of  the  Orr -Sommerfeld  equation,  are  also  of  interest  and 


have  been  studied  by  Gaster  and  Jordinson  in  connection  with  the  evolution  of 


initial  isolated  disturbances  in  the  boundary  layer.  The  neutral  curves, ,on 


which  Im  a = 0,  are  identical  to  the  neutral  curves  of  the  temporal  eigenvalues, 


but  we  do  not  have  such  a detailed  knowledge  of  the  structure  of  a elsewhere 


in  the  complex  plane  as  that  provided  by  Mack  in  the  temporal  case.  We 


can  see  from  Eq.  (4),  however,  that  there  is  a continuous  spectrum  of  eigen- 


2 2 

values  which  are  obtained  by  setting  0^  = -p  and  requiring  p to  be  real. 


The  corresponding  values  of  a are  defined  by 


2 2 
a + iaR  = - p + iwR 


This  spectrum  may  also  be  reinterpreted  in  this  instance  in  terms  of  the 


spectrum  of  the  diffraction  equation.  We  write 


+ = e'1Wt  X(x,  y) 
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and,  when  y ~1,  x » 1,  expand  X in  the  form 

OQ 

la.x—  -s 

X ~ e 2-/Xn(y)x  n (29) 

n = 0 

2 1/2 

where  2a^  = -iR  + i(R  - 4iwR)  , Xn  is  a function  of  y only,  Sn  is  a mono- 
tonic increasing  sequence  of  real  numbers  and  Xq  satisfies 

together  with  the  boundary  conditions  Xq(0)  = X^(0)  = 0.  We  now  forbid  Xq  to 
grow  exponentially  with  y as  y-*oo,  and  then  when  y » 1 

-a,y 

y » D e +Ey  + F (31) 

A-o  o o'  o 

where  D ,E  ,F  are  constants  of  which  E can  be  expected  to  be  non-zero, 
o o o o 

When  y,  x » 1 we  write 


where  t is  monotone  increasing,  X ” ~ a[  r cos  0 = x,  and  r sin  0 = y. 
The  differential  equation  satisfied  by  is 


j_  H __  /V 


(33) 


and  since  G —*•0  as  r— -<x> 
n 


r/R2  \l/2  ■ 

G = (L  cos  m6  + M sin  m6)K  (-3—  - iwR]  r 
o / j m mm  \ 4 / 


where  K is  a Bessel  function  of  the  second  kind.  The  matching  of  Eq.  (29 
m 

as  y—  oo  with  Eq.  (32)  as  y/x— 0 requires  that 

-ay 

X = D e 1 , L =0,  mM  = E (35) 


and  s is  an  integer.  At  this  stage  m remains  indeterminate, 
n 

A parallel  pair  of  expansions  can  be  written  if  x is  large  and  negative 

2 1/2 

except  that  2 must  be  replaced  by  2a ^ = -iR  - i(R  - 4iwR ) and  results 
similar  to  Eq.  (35)  are  recovered.  (Recall  that  negative  x refers  to  stations 
upstream  of  xq  and  not  of  the  leading  edge.  ) It  may  be  expected  that  the 
expansions  for  x » 1 and  x « -1  when  y » 1 are  analytic  continuations  of 
one  another.  Since  the  expressions  for  Gq  in  Eq.  (32)  hold  for  all  G provided 
r » 1,  we  would  then  infer  that  sin  rmr  = 0 also,  in  which  case  m is  an  integer 
and  the  resulting  discrete  spectrum  may  be  interpreted  as  an  integration  over 

the  continuous  spectrum  defined  by  Eq.  (27). 

The  model  equation,  Eq.  (1),  also  has  a continuous  spectrum  of  eigen- 
values, defined  by  Eq.  (27),  which  may  be  reinterpreted  in  terms  of  the 
spectrum  of  the  diffraction  equation.  Further,  it  has  a readily  identifiable 
family  of  discrete  eigenvalues  obtained  by  requiring  R andacm  = w to  be  real 
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and  by  solving  Eq.  (3)  for  a subject  to  the  condition  Re  0^  > 0.  Denoting  a 
typical  member  of  this  family  by  aM-  we  may  easily  show  that  when  R » 1 


? 22  -3/4 

a " = iwR  + n it  + 0(R  ' ) 

M 


where  n is  an  integer.  The  two  solutions  are  presumably  relevant  to  the 
regions  ±x  » 1 and  correspond  to  waves  propagating  upstream  and  down- 
stream. When  R is  finite  the  values  of  aM  must  be  found  numerically.  In 
Figs.  2 and  3 we  display  the  variations  of  aM  with  R for  the  first  two  eigen- 
values; the  wave  velocity  is  positive  in  Fig.  2 and  negative  in  Fig.  3.  These 
eigenvalues  terminate  as  R decreases  when  Re  6^  = 0;  when  the  waves  propa- 
gate downstream  this  occurs  at  R = 6.82  and  14.3,  respectively,  for  the  first 
two  of  them.  For  negative  wave  velocity  the  terminal  values  of  R are 
changed  to  1.  36  and  4.  17. 

Thus  it  follows  that  the  spectrum  of  Eq.  (1)  consists  only  of  a continuous 

family  if  R <1.  36,  which  may  be  interpreted  as  a discrete  spectrum  using 

diffraction  theory.  At  R = 1.  36  a branching  occurs  and  for  R > 1.  36  there  is 

another  eigenvalue  discrete  from  this  family  in  which  Im  <0.  At  R 4.  I i 

another  eigenvalue  branches  off  so  that  for  R > 4.  17  there  are  two  discrete 

eigenvalues,  each  with  Im  < 0,  and  so  on.  We  infer  that  for  general 

R > 6.  82  the  spectrum  of  consists  of  a continuous  family  defined  bv 

Eq.  (27)  together  with  a finite  number  of  discrete  values,  some  of  whi»h  have 

Im  ax  . > 0 and  some  Im  ax,  < 0.  The  number  in  each  set  is  rough  1 y pr  >p*>  r - 
M M 

tional  to  R when  R is  large. 
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Figure  2.  Variation  of  Re  am  and  Im  am  with  R corresponding  to 
waves  propagating  downstream  for  the  first  two  dis- 
crete spatial  eigenvalues  of  the  model  equation  when 
w = 1.  They  branch  off  the  continuous  spectrum  when 
R = 6.  82  and  14.  3,  respectively. 
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d Im  a with  R corresponding  to 
stream  for  the  first  two  discrete 


Figure  3.  Variation  of  Re  am  and  Im  a, 
waves  propagating  ups  I 
spatial  eigenvalues  of  the  model  equation  when  u>  = 1, 
They  branch  off  the  continuous  spectrum  when 
R = 1.  36  and  4.  17,  respectively. 
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5.  THE  DISCRETE  SPATIAL  SPECTRUM  OF  THE 
BLASIUS  BOUNDARY  LAYER 

6 

There  are  at  present  no  calculations  of  the  discrete  spatial  spectrum 

2 

for  the  Blasius  boundary  layer  as  extensive  as  those  given  by  Mack  for  the 
discrete  temporal  spectrum.  Thus,  we  have  developed  a numerical  method 
(which  is  described  in  the  Appendix)  for  computing  the  discrete  spatial  eigen- 
values. Since  Re  a > 0 corresponds  to  disturbances  traveling  downstream 
and  is  of  greater  physical  interest  than  those  in  which  Re  a < 0,  the  numerical 
work  has  been  primarily  concerned  with  the  former.  For  Re  a > 0 we  have 
obtained  the  variation  with  R of  the  two  principal  eigenvalues  corresponding 
to  u)  = 0.  0397,  and  the  results  are  depicted  in  Fig.  4.  The  first  eigenvalue 
springs  from  the  continuous  spectrum  at  R = 10.5;  the  second  at  R = 265. 

The  behavior  of  the  primary  eigenvalue  is  quite  different  from  that  of  the 
model  equation  (Fig.  2);  neither  Rea  nor  Im  a is  monotonic  and  indeed  Im  a 
changes  sign  for  451  < R < 3360.  Furthermore,  it  appears  that  the  primary 
eigenvalue  is  approaching  a limit  as  R— oo;  we  find  this  limit  to  be  a a 0.  14  + 

0.  005i.  The  behavior  of  the  first  discrete  spatial  mode  is  similar  to  that 
2 

observed  by  Mack  for  the  first  discrete  temporal  mode;  the  behavior  for 

4 

R— oo  is  in  agreement  with  Lin's  work.  The  physically  most  important  pro- 
perty of  Eq.  (18)  is  the  possibility  of  instability  over  a finite  range  of  R and 
is  manifested  in  the  unique  behavior  of  the  first  mode.  The  second  mode  in 
Fig.  4 and  presumably  the  higher  modes  have  Rea  and  Im  a monotonic  with 
R and  are  very  similar  to  the  model  (Fig.  2). 
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Figure  4.  Variation  of  Rea  and  Im  a with  R corresponding  to  waves  propagating 
downstream  for  the  first  two  discrete  spatial  eigenvalues  of  the  Orr- 
Sommerfeld  equation  when  w = 0.  0397.  They  branch  off  the  continuous 
SDectrum  when  R = 10.  5 and  265,  respectively. 


The  spatial  Or r -Sommerfeld  problem  has  a continuous  spectrum  with 
Rea  < 0 which  describes  upstream  traveling  waves.  We  have  not  studied  the 
corresponding  discrete  spectrum  in  as  much  detail  for  this  case,  but  have 
found  that  for  w = 0.  0397  a discrete  eigenvalue  springs  from  the  continuous 
spectrum  at  R = 2.  1 with  a = -0.035  - 2.  2i.  Thus,  the  discrete  eigenvalues 
appear  at  a lower  Reynolds  number  than  they  did  for  the  downstream  propa- 
gating waves.  The  upstream  propagating  waves  are  of  course  highly  damped. 
The  quantity  Im  -a  for  the  one  point  given  previously  exceeds  by  an  order  of 
magnitude  the  values  shown  in  Fig.  4.  Equation  (27)  shows  that  the  con- 
tinuous spectrum  is  also  highly  damped;  Im  a = -R  + 0 ( 1 /R  ) when  R » 1 and 
P = 0(1). 
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APPENDIX:  DESCRIPTION  OF  THE  NUMERICAL  SCHEME 

To  obtain  the  results  presented  herein  it  was  necessary  to  solve  the 
Or r -Sommerfeld  equation  numerically  for  selected  discrete  and  continuous 
eigenfunctions  and  eigenvalues.  The  numerical  method  is  similar  to  one  used 
by  Orszag  in  that  the  dependent  variables  are  represented  by  an  appropriate 
series  of  Chebyshev  polynomials . Then  the  discrete  eigenfunctions  and 
eigenvalues  may  be  found  by  application  of  the  QR  algorithm  due  to  Francis  , 
and  the  continuous  solutions  may  be  found  by  solving  a system  of  simultaneous 
linear  equations. 

In  the  present  problem  the  domain  of  y is  semi-infinite;  therefore  a 
convenient  argument  for  the  Chebyshev  polynomials  is  a negative  exponential 
function  rather  than  the  usual  linear  function  of  y.  Thus,  we  have  represented 
a typical  dependent  variable  in  the  form 

N 

*(y)  = (e~y/yf)  (37) 

j=o  J J 

where  T is  the  Chebyshev  polynomial  defined  on  the  interval  zero  to  one 

and  y is  a scale  factor.  (The  properties  of  Chebyshev  polynomials  are  given, 

9 

for  example,  by  Fox  and  Parker  . ) 

The  first  step  in  the  numerical  procedure  is  to  solve  the  Blasius  equa- 
tion so  it  may  be  represented  in  the  form  of  Eq.  (37).  This  is  achieved  by 
iterating  the  Weyl  integral  formulation  of  the  Blasius  equation,  as  suggested 
by  Jones  and  Watson.^  With  a 17-term  expansion  (N  = 16)  and  y^  = 4.  Z5  the 
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solution  for  U has  a maximum  error  of  order  10  on  the  interval  0 < y £ oo. 
Since  N is  always  greater  than  or  equal  to  16  for  the  results  presented  herein, 
any  errors  in  the  Blasius  solutions  are  negligible. 

The  Or  r - Somme  rfe  Id  equation  is  solved  by  formulating  the  temporal 
eigenvalue  problem  with  the  a/s  in  Eq.  (37)  as  unknowns.  This  results  in 
a linear  eigenvalue  problem  and  the  QR  algorithm  is  used  to  solve  it.  (The 
method  is  similar  to  the  one  used  by  Orszag'  to  solve  the  Or r -Somme rfeld 
equation  for  plane  Poiseuille  flow.  ) 

For  given  values  of  a and  R,  the  eigenvalues  of  the  Or r -Somme rfeld 
equation  consist  of  a continuous  spectrum  and  a finite  number  of  discrete 
values.  The  numerical  model  has  N-2  eigenvalues,  some  of  which  are  spurious 
and  some  of  which  approximate  the  actual  spectrum.  We  make  two  tests  for 
the  relevance  of  the  solutions.  First  the  value  of  the  eigenvalue  should 

approach  a limit  as  N is  increased;  second  the  series  Eq.  (37)  describing  the 

_ ^ 

eigenfunction  must  be  convergent  and  accurate  to  order  10  or  less.  [That 

is  the  a's  as  j—  N must  be  0(10  ).  ] About  a third  of  the  eigenvalues  can  be 

excluded  on  the  basis  of  the  first  test  since  they  do  not  approach  a limit. 

Over  half  the  eigenvalues  are  very  near  the  continuous  spectrum  and  have 

-4  - 1 

Re  c = 1 - 6 with  c varying  between  10  and  10  in  a typical  case.  The 
quantity  Im  c does  not  approach  a limit  with  N nor  does  the  stream  function 
series  Eq.  (37)  converge.  We  interpret  these  solutions  as  approximate 
representations  of  the  continuous  spectrum  by  the  discrete  numerical  model 
but  recognize  that  this  representation  can  never  be  accurate  because  the 


function  set  is  countable  and  decays  exponentially  while  the  continuous  spec- 
trum is  sinusoidal  at  infinity.  The  remaining  solutions  represent  the  dis- 
crete eigensolutions  of  the  Or r -Somme rfeld  equation.  In  the  Reynolds  number 
range  in  which  instability  occurs,  rather  modest  values  of  N will  result  in 

three-figure  accuracy  for  the  primary  discrete  eigensolution;  both  larger 

2 5 

and  smaller  values  of  R require  an  increased  N.  For  example,  with  R = 10  , 

a - 0.  115,  N = 16,  and  y^  = 4.25,  the  maximum  error  in  the  stream  function 
_3 

is  0(10  ) compared  to  unity  on  the  semi-infinite  interval.  To  obtain  the 

curves  of  the  first  two  spatial  modes  shown  in  Fig.  4 we  have  used  values  of 

N up  to  50;  the  value  of  yg  has  also  been  varied  to  minimize  N.  (The  optimum 

value  of  v varies  directly  with  some  characteristic  thickness  of  the  eigen- 
’ e 

solutions  and  therefore  decreases  with  R.  ) 

Solutions  to  the  temporal  eigenvalue  problem  are  easily  obtained:  R 
and  a real  value  of  a are  specified,  and  in  some  cases  more  than  one  mode 
may  be  obtained  from  one  calculation.  The  spatial  eigenvalue  problem  is 
more  complicated:  R is  specified  and  we  iterate  the  complex  value  of  a 
until  u>  has  a specified  real  value.  This  procedure  works  well,  converging  in 
two  or  three  iterations.  However,  the  first  one  or  two  points  of  the  higher 
modes  are  more  difficult  to  obtain  because  there  is  no  convenient  starting 
point  as  there  is  with  the  first  mode  which  has  a and  u>  both  real  at  the  stability 
boundary. 

The  representation  of  the  eigenfunctions  as  polynomials  in  negative 
exponentials  is  well  suited  to  functions  of  the  boundary-layer  type.  However, 
as  the  discrete  eigensolutions  approach  the  intersection  point  with  the 


continuous  spectrum,  the  thickness  increases  and  the  numerical  scheme 
described  herein  fails.  In  order  to  locate  these  points  a numerical  method 
for  generating  solutions  in  the  continuous  spectrum  is  used.  Approximate 
solutions  of  the  Or r -Somme rfeld  equation  in  the  continuous  spectrum  have 
been  obtained  by  Rogler  and  Reshotko  ; however,  they  considered  only  the 
case  p = a[cf.  Eq.  (21)]  and  they  neglected  the  imaginary  part  of  c.  The 
solution  technique  used  herein  is  similar  to  that  of  Rogler  and  Reshotko; 
that  is,  we  subtract  the  oscillatory  part  of  the  solution  from  the  dependent 
variable  and  thereby  obtain  an  inhomogeneous  version  of  the  Or  r -Sommerfeld 
equa  tion. 

®(y)  = 'l'(y)  - C sin  py  - D cos  py  (38) 

[in  some  cases  we  found  that  accuracy  was  increased  by  also  subtracting 
E exp  (-ay)  from  ^ in  Eq.  (38).]  The  solution  procedure  is  as  follows: 
Substitute  Eq.  (38)  into  the  Or r -Somme rfeld  equation  and  the  boundary  con- 
ditions and  specify  R,  p,  and  a (or  w).  Then  w (or  a)  is  given  by 

o'  + p^  = iR(w-a)  (39) 

and  ® and  the  constants  C and  D are  unknowns  defined  by  a linear  system  of 
equations  (not  an  eigenvalue  problem).  We  solve  this  system  using  a 
Chebyshev  expansion  of  the  form  Eq.  (37),  but  that  is  not  essential. 


The  intersection  points  are  those  points  at  which 


|C|  = j D j ; arg  D - arg  C = tt/2  (40) 

We  locate  these  points  by  generating  the  curve  as  a function  of  R and  p upon 
which  the  first  condition  of  Eq.  (40)  is  satisfied.  The  point  on  the  curve 
where  both  of  Eqs.  (40)  are  satisfied  is  an  intersection  point. 
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